SPARTA (7 Jul 2014)
# advect particles via VSS collisional flow on a uniform grid
# particles reflect off global box boundaries

variable            x index 10
variable            y index 10
variable            z index 10

variable            lx equal $x*1.0e-5
variable            lx equal 100*1.0e-5
variable            ly equal $y*1.0e-5
variable            ly equal 100*1.0e-5
variable            lz equal $z*1.0e-5
variable            lz equal 100*1.0e-5

variable            n equal 10*$x*$y*$z
variable            n equal 10*100*$y*$z
variable            n equal 10*100*100*$z
variable            n equal 10*100*100*100

seed	    	    12345
dimension   	    3
global              gridcut 1.0e-5

boundary	    rr rr rr

create_box  	    0 ${lx} 0 ${ly} 0 ${lz}
create_box  	    0 0.001 0 ${ly} 0 ${lz}
create_box  	    0 0.001 0 0.001 0 ${lz}
create_box  	    0 0.001 0 0.001 0 0.001
Created orthogonal box = (0 0 0) to (0.001 0.001 0.001)
create_grid 	    $x $y $z
create_grid 	    100 $y $z
create_grid 	    100 100 $z
create_grid 	    100 100 100
Created 1000000 child grid cells
  parent cells = 1
  CPU time = 0.171848 secs
  create/ghost percent = 79.8856 20.1144

balance_grid        rcb part
Balance grid migrated 0 cells
  CPU time = 3.59247 secs
  reassign/sort/migrate/ghost percent = 3.03936 0.234545 0.592691 96.1334

species		    ar.species Ar
mixture		    air Ar vstream 0.0 0.0 0.0 temp 273.15

global              nrho 7.07043E22
global              fnum 7.07043E6

collide		    vss air ar.vss

create_particles    air n $n
create_particles    air n 10000000
Created 10000000 particles
  CPU time = 2.27579 secs

stats		    10
compute             temp temp
stats_style	    step cpu np nattempt ncoll c_temp

# first equilibrate with large timestep to unsort particles
# then benchmark with normal timestep

timestep 	    7.00E-8
run                 30
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 1029.97 1029.97 1029.97
  grid      (ave,min,max) = 175.552 175.552 175.552
  surf      (ave,min,max) = 0 0 0
  total     (ave,min,max) = 1205.52 1205.52 1205.52
Step CPU Np Natt Ncoll temp 
       0            0 10000000        0        0    273.15661 
      10    63.774604 10000000  9632983  7045171    273.15661 
      20    140.59544 10000000  9895865  7057044    273.15661 
      30    223.07336 10000000 10031028  7061938    273.15661 
Loop time of 223.073 on 1 procs for 30 steps with 10000000 particles

Particle moves    = 300000000 (300M)
Cells touched     = 1486854861 (1.49B)
Particle comms    = 0 (0K)
Boundary collides = 11989295 (12M)
Boundary exits    = 0 (0K)
SurfColl checks   = 0 (0K)
SurfColl occurs   = 0 (0K)
Collide attempts  = 290489215 (290M)
Collide occurs    = 210911693 (211M)
Particles stuck   = 0

Particle-moves/CPUsec/proc: 1.34485e+06
Particle-moves/step: 1e+07
Cell-touches/particle/step: 4.95618
Particle comm iterations/step: 1
Particle fraction communicated: 0
Particle fraction colliding with boundary: 0.0399643
Particle fraction exiting boundary: 0
Surface-checks/particle/step: 0
Surface-collisions/particle/step: 0
Collision-attempts/particle/step: 0.968297
Collisions/particle/step: 0.703039

Move  time (%) = 125.983 (56.476)
Coll  time (%) = 87.2 (39.0903)
Sort  time (%) = 9.58561 (4.29707)
Comm  time (%) = 0.00011611 (5.20501e-05)
Outpt time (%) = 0.304642 (0.136566)
Other time (%) = 6.81877e-05 (3.05674e-05)

Particles: 1e+07 ave 1e+07 max 1e+07 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Cells:      1e+06 ave 1e+06 max 1e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0
GhostCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
EmptyCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
timestep 	    7.00E-9
run 		    100
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 1029.97 1029.97 1029.97
  grid      (ave,min,max) = 175.552 175.552 175.552
  surf      (ave,min,max) = 0 0 0
  total     (ave,min,max) = 1205.52 1205.52 1205.52
Step CPU Np Natt Ncoll temp 
      30            0 10000000 10031028  7061938    273.15661 
      40    23.372189 10000000   902713   704702    273.15661 
      50    47.152074 10000000   917836   705180    273.15661 
      60    71.041471 10000000   930209   705855    273.15661 
      70    94.995363 10000000   939663   706542    273.15661 
      80    118.98444 10000000   947092   705254    273.15661 
      90    143.10531 10000000   953647   706201    273.15661 
     100    167.31628 10000000   959994   705970    273.15661 
     110    191.57303 10000000   964369   705586    273.15661 
     120    215.86028 10000000   969484   706558    273.15661 
     130    240.22351 10000000   971927   705583    273.15661 
Loop time of 240.224 on 1 procs for 100 steps with 10000000 particles

Particle moves    = 1000000000 (1B)
Cells touched     = 1395629264 (1.4B)
Particle comms    = 0 (0K)
Boundary collides = 3994827 (3.99M)
Boundary exits    = 0 (0K)
SurfColl checks   = 0 (0K)
SurfColl occurs   = 0 (0K)
Collide attempts  = 93687763 (93.7M)
Collide occurs    = 70175527 (70.2M)
Particles stuck   = 0

Particle-moves/CPUsec/proc: 4.16279e+06
Particle-moves/step: 1e+07
Cell-touches/particle/step: 1.39563
Particle comm iterations/step: 1
Particle fraction communicated: 0
Particle fraction colliding with boundary: 0.00399483
Particle fraction exiting boundary: 0
Surface-checks/particle/step: 0
Surface-collisions/particle/step: 0
Collision-attempts/particle/step: 0.0936878
Collisions/particle/step: 0.0701755

Move  time (%) = 110.289 (45.911)
Coll  time (%) = 97.0061 (40.3816)
Sort  time (%) = 31.9155 (13.2858)
Comm  time (%) = 0.000422955 (0.000176067)
Outpt time (%) = 1.01215 (0.421338)
Other time (%) = 0.000219822 (9.15072e-05)

Particles: 1e+07 ave 1e+07 max 1e+07 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Cells:      1e+06 ave 1e+06 max 1e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0
GhostCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
EmptyCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0

